Supersymmetric isospectral formalism for the calculation of near-zero energy states: 
application to the very weakly bound ^He trimer excited state 
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We propose a novel mathematical approach for the calculation of near-zero energy states by solv- 
ing potentials which are isospectral with the original one. For any potential, families of strictly 
isospectral potentials (with very different shape) having desirable and adjustable features are gen- 
erated by supersymmetric isospectral formalism. The near-zero energy Efimov state in the original 
potential is effectively trapped in the deep well of the isospectral family and facilitates more accurate 
calculation of the Efimov state. Application to the first excited state in *E[e trimer is presented. 

PACS numbers: 03.75.Hh, 31.15.Ja, 03.65.Gc, 03.75.Nt 



I. INTRODUCTION 



It was proposed by Efimov in 1970 that if two spinless 
neutral bosons interact resonantly then the addition of 
a third identical particle leads to the appearance of an 
infinite number of bound three-body energy levels [l|, Q . 
This occurs simultaneously with the divergence of the 
s-wave scattering length (os). associated with appear- 
ance of an additional zero-energy two-body bound state. 
Hence highly exotic Efimov states appear when there 
is a zero or near-zero energy two-body bound state. 
For a long time there was no clear signature of Efimov 
states in any naturally occuring trimer system. Efimov 
states are not possible in atomic systems due to the 
long range Coulomb interaction, however it may exist 
in the system of spinless neutral atoms. Even though 
the Efimov effect was predicted four decades ago (T], , 
evidence of its existence in ultracold caesium and potas- 
sium trimers has been experimentally established only 
very recently Q. However, these trimers are obtained 
by manipulating two-body forces through Feshbach 
resonances and are not naturally occuring. Therefore, 
it is of great interest to search for the Efimov effect 
in a naturally occuring trimer, like *He trimer. So far 
no experimental confirmation has been reported. The 
near-zero energy {Eq ^ 1 mK) bound state (which is also 
the only bound state) of '^He dimer opens the possibility 
of the existence of an Efimov-like state in "'He trimer. 
Several authors remarked that the ^He trimer may 
be the most promising candidate. Earlier theoretical 
calculations show that the trimer has a L =0 ground 
state at 126 mK and a near-zero energy excited state 
(- 2mK) [M3. The excited state has been claimed to 
be an Efimov state. A controversy arises from the fact 
that the number of Efimov states is highly sensitive to 
the binding energy of the dimer and even a very small 
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decrease of the strength of two-body interaction makes 
the system unbound. Strikingly, it also disappears when 
the two-body interaction strength is increased. However 
in contrast with theoretical investigations, no evidence 
of Efimov trimer has been found experimentally [l5| . 
In the experiments, ^He trimer has been observed in 
its ground state only. No experimental evidence of the 
excited state has been reported so far. 

In principle "^He trimer may be considered as a very 
simple three-body system consisting of three identical 
bosons. But its theoretical treatment is quite difficult. 
First, the He-dimer potential is not uniquely known. 
Very small uncertainities in the dimer potential may lead 
to different conclusions. Secondly, the strong short-range 
interatomic repulsion in the He-He interaction causes 
large numerical errors. As *He systems are strongly 
correlated due to large ''He-^He repulsion at short 
separation, the effect of interatomic correlation must be 
taken properly into account. 

In the present communication, we revisit the problem 
using a correlated basis function known as potential 
harmonics (PH) basis which takes care of two-body 
correlations In order to include the effect of 

highly repulsive He-He core, we multiply the PH basis 
with a suitable short-range correlation function which 
reproduces the correct short-range behavior of the 
dimer wavefunction. Although this correlated PH 
basis (CPH basis) correctly reproduces the dimer and 
trimer properties, we could not find any Efimov like 
state in trimer with the actual dimer interaction jl7| . 
We point out that the calculation of such a near-zero 
energy excited state in the shallow and extended trimer 
potential may involve severe numerical errors and we 
may miss it. Thus an alternative accurate procedure is 
desirable. Here, we apply the supersymmetric isospectral 
formalism for an accurate treatment. For any given 
potential, families of strictly isospectral potentials, with 
very different shape but having desirable and adjustable 
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Fig. 1 (Color online) The eilective potential ujo{r) (red solid 
curve) and isospectral potentials d)o(j", A) corresponding to 
two values of A: A = 0.00005 (green dashed curve) and 
A = 0.00002 (blue dotted curve) for the ''He trimer. All 
energies are in mK and r in a.u. The horizontal line 
indicates the energy of the first excited state in cjoir)- 



features ar e g enerated by supersymmetric isospectral 
formalism [18[. The near-zero energy bound state will 
be more effectively bound in the deep narrow well of 
the isospectral potential and will facilitate an easier and 
more accurate calculation of the near-zero energy excited 
state. Following the steps of supersymmetric quantum 
mechanics [isj . for any given potential wo(f), one can 
construct a class of potentials ujQ{r,{X}), where {A} 
represents a set of one or more continuously variable real 
parameters. The potential wq is isospectral in the sense 
that ujQ and ujq have identical spectrum, reflection and 
tranmission coefficients. For simplicity wc consider only 
one parameter (A) family of isospectral potentials. Wc 
will see later that A can take real values — oo < A < — 1 
and < A < cxj. For A — ?> oo, one gets back the original 
potential. Although the set of isospectral potentials are 
strictly isospectral with the original potential, they have 
different shapes depending on the parameter A [1^. In 
Fig. 1, we demonstrate how an original potential, ujQ{r) 
shown by the solid (red) curve, changes in the isospectral 
potential u}Q{r, X) for two values of the parameter A, 
viz. A = 0.00005 (green dashed curve) and A = 0.00002 
(blue dotted curve). We introduce this figure here 
for a qualitative understanding of the features of the 
isospectral potentials. A complete discussion of how such 
isospectral potentials are calculated will be presented 
in Sections IIC and IIL Although all three potentials 
produce identical energy spectrum, their shapes are seen 



to be very different. The original potential [wo('')] has a 
shallow and wide well with a short range repulsion. By 
contrast, both the isospectral potentials have a deep and 
narrow attractive well (NAW) at smaller r, while the 
long range part does not differ much from uJo{r). One 
can also notice that as A decreases, the narrow attractive 
well becomes deeper, while the intermediate barrier 
becomes higher. Hence, the near-zero energy (indicated 
in Fig. 1 by a horizontal line) excited state in uio{r) will 
be at the same energy value in the isospectral potentials, 
but now that state will lie deep within the NAW. Thus, 
while this state is weakly bound and spatially extended 
in the original potential, it will be strongly bound and 
well localized in the NAW of the isospectral potential, 
at the same energy. Clearly, computation of the wave 
function and its energy (equal to the energy of the first 
excited state in the original potential) will be easier. 
Furthermore, this state becomes more strongly bound, 
as A decreases. In general, as A decreases continuously 
from oo to 0, wq starts developing a local minimum 
which shifts towards r=0 and becomes deeper and 
narrower. Consequently, a shallow potential transforms 
into one having a narrow and deep potential well near 
the origin in the isospectral potential. The surface 
barrier also becomes high. Such interesting properties of 
Wo can be useful to solve near-zero energy states. Such 
a state lies near the top of the original potential well 
and its wave function is spatially very extended, while 
the isospectral state lies well within the NAW. Hence 
the latter is strongly bound and well localized within 
the narrow well of Vi- The parameter A controls these 
features and a suitable optimum value can be chosen 
(see later). 

Thus our approach consists of two steps. First, to ap- 
ply a correlated quantum many-body theory for a highly 
correlated system like *He-trimer and second, to use the 
isospectral formalism for the accurate determination of 
the first excited state, whose energy is just a few mK. 

The paper is organized as follows. In Section II, 
we present a brief review of the correlated potential 
harmonics expansion method, choice of potential and 
the isospectral formalism. Section HI presents the 
results of our numerical calculation. Finally we draw 
our conclusions in Section IV. 



II. THEORETICAL PROCEDURE 

A. Correlated Potential harmonics expansion 
method 

The Hamiltonian for a system of {N + 1) atoms ( each 
of mass m) and interacting via two-body potential has 
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the form 
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(1) 



where V{xi ~ xj) = F(rij) is the He- He two-body poten- 
tial described later. The relative motion in the standard 
Jacobi coordinates, {(fi, . . . , Ctv}, is given by [ll] 



— J2% + yiCi,---XN)^E 



*(Ci,...,Civ) = 0. 

(2) 



We decompose ^ in Faddeev components 



TV -1-1 

E 
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(3) 



where -tpij (x) is the two-body Faddeev component for the 
(ij) partition. In potential harmonics expansion method 
(PHEM), we expand (ij)-Faddeev component in the com- 
plete set of potential harmonics, {'p2'ff-(-i(^JV ')}' appro- 
priate for (ij)-partition [TgI. Hoj 



r). 



(4) 



K 



PH basis is the subset of hyperspherical harmonics (HH), 
which is sufficient for the expansion of the two-body po- 
tential V{fij), for the particular {ij) partition [l^. This 
PH function depends only on and a global length r 
(called hyperradius) . Hence in Eq. ([3]) includes only 

(ii) 

long range two-body correlations, fijv represents the 
full set of hyperangles for the (ij) partition. The po- 
tential harmonics basis function for the (ij)-partition is 
given by [16] 



P2K+1 iS^ij ) 



,„K)Wp^^,^,(</.)yo(i?-3), (5) 



where 



= cos ^{rij/r), Yim{i^ij) is a spherical har- 
monic, ujij being the polar angles of Vij, P2k+i('^) 
is expressed in terms of Jacobi polynomials and j^o is 
the lowest order hyperspherical harmonic in {^N — 3) 
dimensional hyperangular space (hence a constant) [loj . 
The basic assumption in Eq. (3) is that only two-body 
correlations are important and higher-body correlations 
are disregarded. Thus in the (ij)-Faddeev component, 
where only the (ij)-pair interacts, tpij is independent of 
the coordinates of all particles except . Hence we can 
freeze the contributions coming from (N — 1) remaining 
spectators. Thus the contribution to the orbital angular 
momentum and the grand orbital quantum number 
comes only from the interacting pair and the 3A^ di- 
mensional Schrodinger equation reduces effectively to a 
four dimensional equation. The relevant set of quantum 
numbers, associated with the hyperangles, are three. 
These are orbital I, azimuthal m and grand orbital 
2K + I for any N. 



So far we have disregarded the effect of strong short 
range correlation in the PH basis. The He-He potential 
becomes suddenly very strongly repulsive below a certain 
value of interatomic separation. This causes a very strong 
short range two-body correlation in the many body wave 
function. We introduce this correlation function in the 
expansion basis and call it as CPH basis [l^] . 



2K+1 



correlated 



= yz™(c.,,) wp^^^,(0) 

3^o(3iV-3)?7(r,,), 



(6) 



where rj{rij) is the short range correlation function. It 
is practically zero for very small Vij (rij < Vc) where Tc 
is the size of the repulsive core. Its role is to enhance 
the speed of convergence of the expansion basis. We ob- 
tain r]{rij) as the zero energy solution of (ij)-pair relative 
motion in the potential V{rij), 



1 d 



dri{r 



dr. 



+ V{nMnj) = 0. (7) 



Replacing PH by CPH in Eq. (4) and taking projection of 
the Schrodinger equation on the PH basis of a particular 
partition, a set of coupled differential equation (CDE) is 
obtained H 



- — TT + ^{>C(£ + l)+4K(/f + a + /3 + l)} 



- E\uKi{r) 

+ fKlVKK'{r)fK'lUK'l{r 



(8) 



0, 



K' 



where C = I 



3N-3 



3W-5 



^— , ^ — , P = I + ^ and fxi is 

a constant representing the overlap of a PH correspond- 
ing to a particular partition with the sum of PHs of all 
partitions [20| . K is the hyperangular momentum quan- 
tum number. The correlated potential matrix element 
Vkk' {t) is now given by (20l | 



^Pu^P^~i / Pf{z)v{r 



l + z 



(9) 



Here h'^ and Wi{z) are respectively the norm and 
weight function of the Jacobi polynomial P^^ [l6| . 



B. Choice of He-He potential 

For a realistic calculation, one needs an accurate He- 
He interaction potential. Several sophisticated He- He po- 
tentials have been proposed [2l|. Among these, the com- 
monly used ones are: Tang, Tonnies and Yiu (TTY) (2^ . 
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LM2M2 [23|, and HFD-HE2 [24| potentials. These poten- 
tials reproduce all known two-body He-He data. In the 
present work, we select the more popular and sophisti- 
cated TTY potential. This potential has the form [lTI.[22| 



V{x) = A[Ve,{x) + Vmsp{x)1 



(10) 



where x represents the intcrparticle distance. The 
part Vex has the form Vex{x) = DxPe~^'^^ wih p = 



^ — 1. The other part Vdisp is given as Vdisp{x) 

- -EiL3^2„/2„(x)x-2". 



The coefficients C271 
are calculated using the recurrence relation C2n = 

(§^)'c'2„-6; C'e = 1.461, Cs = 14.11, Cw = 183.5, 

A = 315766.2067 K, D= 7.449 and 7 = 1.3443 {a.u.y^. 

The function /2„ is given by f2n{x) = 1-e''"' J2k=o TT- 
with b{x)= 27 - f. 



For our numerical solution, the set of CDEs [Eq. (8)] 
is solved by hyperspherical adiabatic approximation 
(HA A) [2^. In HA A, one assumes that the hyperra- 
dial motion is slow compared to the hyperangular mo- 
tion. The effective potential for the hyperradial motion 
(obtained by diagonalizing the potential matrix together 
with the diagonal hypercentrifugal repulsion for a fixed 
value of r) is obtained as a parametric function of r. 
We choose the lowest eigenpotential {u!o{r)) as the effec- 
tive potential. Thus in HAA, energy and wavefunction 
are obtained approximately by solving a single uncoupled 
differential equation 



, +ujQ{r)-E 



Coir) = 0, 



(11) 



subject to appropriate boundary conditions on Co('')- 
The principal advantage of the present method is 
two-fold. Firstly, the CPH basis set correctly takes 
care of the effect of strong short range correlation 
produced by the He-He interaction. Secondly, the use 
of HAA basically reduces the multidimensional problem 
to an effective one dimensional problem introducing the 
effective potential. The effective potential {ujo{r)) gives 
a clear qualitative as well as quantitative picture. For 
our numerical calculation we investigate the I = state 
and truncate the CPH basis to a maximum value K = 
Kmax, requiring proper convergence. 

The ground state properties of ^He dimer is obtained as 
a numerical solution of two-body Schrodinger equation by 
Runga-Kutta algorithm. The dimer energy using TTY 
potential, as well as the results from other references are 
presented in Table I. 



Although in our earlier calculation [l^ we reported the 
dimer and trimer ground state properties, we include 
them here for completeness of the discussion. Calculated 
rms value of is 98.596 a.u. The extremely small 



TABLE I: The ''He-dimer energy using TTY potential. 



Expt. 


ed(mK) 


1 1+0.3 


present method 


-1.254 


mK 


DMC [23] 


-1.243 


[27] 


Other [11] 


-1.309 




Other [8] 


-1.313 



binding energy of the dimer and the large spatial 
extension of the ground state wave function imply that 
the ground state of helium dimer is very loosely bound. 
The trimer ground state energy, as well as the results 
obtained in earlier investigations by other authors are 
presented in Table II for different potentials. The r.m.s. 
value of hyperradius is 21.389 a.u. Thus our correlated 
PH basis successfully reproduces the energy values which 
are in very close a gree ment with other sophisticated 
calculations d-H [Ullillil . 



TABLE 11: The absolute value of *He-trimer ground state 
energy (in units of mK) obtained by different methods. 



Potential 


PHEM 


variational 


Faddeev 


Adiabatic 


DMC 


TTY 


125.51 


126.40[24] 


126.4 [7] 




125.46 [25] 


LM2M2 


126.37 


126.40 [25] 


126.4 [7] 


125.2 [6] 




HFDHE2 


120.28 




117.1 [7] 


98.1 [5] 





Although we have done detailed calculation of ^He 
trimer ground state energy, we failed to obtain any trimer 
excited state with this choice of two-body interaction. 
Very recently we have analyzed in details the behavior of 
*He trimer excited states as a function of pairwise inter- 
action 



VHe- 



He 



6Vn 



TTY 



(12) 



where 6 controls the strength of two-body interaction. 
(5 =1 is the physical value of the dimer interaction and 
(5=0 corresponds to free particle limit when neither 
two-body nor three-body bound states appear. We 
found that by increasing S from a small value, the trimer 
starts to support an excited state at 6 =0.978 (T7| . 
Binding energy of the excited state gradually increases 
with increase in (5, attains its maximum value at S = 
0.984, then it decreases gradually and disappears which 
indicates dissociation to trimer fragments as a dimer and 
monomer. Thus the disappearance of the first excited 
state due to both increasing and decreasing S clearly 
show that the state is an Efimov state. But the value 
of S for which this happens, does not correspond to the 
actual physical dimer interaction (S = 1). The Efimov 
property of other excited states are also discussed in our 
earlier study |17| . However we did not claim that the 
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trimer excited state does not exist at all, as there is a 
large body of work in this direction |4|-[l3{. So this is 
still an ongoing issue as there is considerable controversy 
in the earlier discussions found in Refs Thus our 

earlier work could not resolve the problem. It may be 
the limitation of our basis set which is unable to produce 
such an elusive state. So it needs further study. 



C. Isospectral formalism 

For an accurate determination of the elusive near- 
zero energy excited state, we use the supersymmetric 
isospectral formalism mentioned earlier. Note that this 
state lies near the top of the potential well uJo{r) and 
is very extended spatially. Consequently, a convergent 
calculation requires a very large number of hyperspher- 
ical partial waves (corresponding to very large Kmax)- 
On the other hand, the trimer ground state lies near 
the bottom of this well, is well bound and localized; it 
therefore converges relatively easily. Now the isospectral 
potential u)i:)[\^r) with sufficiently small positive A will 
have a deep well near the origin, followed by a high 
barrier (see below). The near-zero energy state sought 
after is now strongly trapped within this narrow well, i.e. 
it is sharply localized. Hence its convergence is achieved 
easily enough. The isospectral potential is obtained in 
terms of the trimer ground state wave function (see 
below) , which by the above argument is determined with 
a fair degree of confidence and it offers easier calculation 
of the excited state. 

Since isospectral formalism is not a common topic, we 
briefly explain how an isospectral potential is obtained 
for a given potential [l^. For the given potential wo(r) 
having a normalized ground state Co(^) with energy i?o, 
[see Eq. (11)], a superpotential W{r) is defined as 



V™Co('') 
Then it can easily be seen that 



(13) 



satisfied by the unknown function W{r) (called Riccati 
equation). With this W{r), we can get back wo(?'), using 
Eq. (14). But solution of the non-linear Eq. (15) is not 
unique. For simplicity, we use the units in which ^ — \. 
Then the most general solution is [l8| 



W{r,\) = W{r) 



dr 



ln|/o(r)+A|, 



(16) 



where A is an integration constant and W{r, A) is a func- 
tion of r. parametrically dependent on A. /o('') is given 

by 



(17) 



J 

Then the family of potentials ujo{r^ A) given by 



wo(?', A) - -Bo 
a)o(r, A) 



W'^{r, A) - W'{r, A) 



'' dr- 



■ln|/o(r) + A| 



(18) 



for all allowed values of A (see below), has the same 
partner u}^'^\r). Hence the family of potentials given 
by Eq. ([T8| aU have identical spectrum. Since Co('') is 
normalized, < /o('') < 1- Hence from Eq. (|16|) . one 
notices that the interval — 1 < A < is not allowed. 
For all other values of A, iliQ{r^\) is strictly isospectral 
with ujQ{r) [ll, [s^l- For A = oo, uji^{r,\) becomes a;o(r). 
For small positive values of A, the isospectral potential 
Liji^{r,\) develops an attractive well near the origin. As 
A — >■ 0, this well becomes deeper followed by a high 
barrier, before the shallower part. Note that the ground 
state of wo(r), viz. Co(''), is necessary to calculate the 
family of isospectral potentials. In the present study the 
ground state of oJo{r) is fairly accurately calculated. We 
will sec in the next section that different eigen states are 
transformed differently. Therefore this transformation 
cannot be considered as a generalized rotation in the 
Hilbert space. Thus the transformation from the original 
Hamiltonian to the isospectral Hamiltonian is different 
from a standard unitary transformation, even though 
the entire energy spectrum is preserved. 



uja{r) -Ea^ W^{r) - ^W'ir 



(14) 



and one can define a supersymmetric partner potential 
w(^)(r) through 

uj^^\r)^Eo = WHr) + ^W'{r), (15) 



such that uJo{r) and its partner have the same energy 
spectra, except that the ground state of u!o{r) is absent 
in the spectrum of its partner (Tsj . 

Now, for the given partner potential a;^^^(r), Eq. (15) 
can be considered as a non-linear differential equation 



III. RESULTS AND DISCUSSIONS 

First we calculate the lowest effective potential u!o{r) 
in the hyperradial space by hyperspherical adiabatic 
approximation (HAA) as stated in Sec. IIB. To calculate 
the energy E and wave function Co('') of the original 
potential we solve the single uncoupled equation, 
Eq. (11), with appropriate boundary conditions. As 
stated earlier, although our calculated ground state 
energy Eq and wave function are in good agreement 
with other calculations, we fail to get the first excited 
state in such a shallow potential. Here the isospectral 
formalism will be an effective technique to calculate 
several isospectral potentials of gradually varying shape, 
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which will facilitate easier calculation of the very weakly 
bound state. We calculate /o('') from Eq. (17) and 
then the isospectral potential for a specific value of A is 
calculated from Eq. (18). 

In Fig. 1, the effective potential uJo{r) calculated by the 
CPH method is shown as a continuous (red) curve. This 
has a a soft repulsion at smaller r, followed by a shal- 
low minimum which supports the ground state at 125.51 
mK. As stated earlier, we fail to get the first excited state 
in this effective potential, by our numerical calculation. 
Next we calculate the isospectral potential ujQ{r,X), for 
chosen values of A. We checked that for large values of A, 
calculated u!o(r,X) is practically indistinguishable from 
u!o{r)- For small values of A, the isospectral potential de- 
velops a deep and narrow well near the origin (left side 
well, LSW) with a barrier (intermediate barrier, IB) sep- 
arating it from the shallow well (right side well, RSW). 
As A decreases towards zero, LSW becomes deeper, nar- 
rower and closer to the origin and at the same time, IB 
becomes higher. 



TABLE III: Parameters of the original and the isospectral po- 
tentials. A = CO corresponds to the original potential. LSW, 
IB and RSW respectively represent the left side well, the in- 
termediate barrier and the right side well. Position values are 
in a.u. and energies are in mK. 



A 


LSW 


IB 


RSW 


left 


left 
min 




^max 


right 
min 


riqht 
min 


oo 










8.999 


-3.238 


0.00010 


4.607 


-1.897 


5.327 


8.333 


9.087 


-2.848 


0.00005 


4.527 


-8.119 


5.137 


10.23 


9.087 


-2.847 


0.00002 


4.457 


-22.133 


4.943 


13.874 


9.085 


-2.846 


0.00001 


4.417 


-40.254 


4.817 


18.019 


9.085 


-2.846 



In Table III, we present numerical values of the parame- 
ters of the original (corresponding to A = oo) and isospec- 
tral potentials wo(r,A) for A = 0.00010,0.00005,0.00002 
and 0.00001. For each isospectral potential, we give 
values of the position {r^i^) and value (w^{^) of the 
minimum of LSW, position {r„^ax) a-nd value {oJmax) 
of the maximun of IB, and position (j'^f^*) and value 

i'^min ) of the minimum of RSW. We have checked by 
numerical calculation that the energy of both the ground 
and the first excited state remains independent of the 
choice of A, as theory predicts. However, very small val- 
ues of A arc not convenient for numerical calculation of 
the energy of the excited state, as numerical errors creep 
in, due to extreme narrowness of the LSW, in which the 
excited state resides. Thus a judicious choice of the value 
of A is necessary. By careful investigations for different 
values of A, we find that the range 0.00002 - 0.00005 for 
the value of A is optimum for minimizing errors. As two 
typical cases demostrating the behavior of isospectral 



potentials, we plot ujQ{r, A) for A = 0.00002 (blue, dotted 
curve) and A = 0.00005 (green, dashed curve) in Fig. 1. 
The very weakly bound first excited state in the original 
shallow potential ijJo{r) now becomes strongly bound 
within the deep and narrow LSW of LUo{r, A). The deep 
well and the adjacent high barrier strongly localizes this 
state in the isospectral potential. We solve Eq. (11), 
with LOoif) replaced by LOoir^X), subject to appropriate 
boundary conditions to calculate the energy of the first 
excited state. Following this isospectral technique, we do 
indeed get a bound first excited state. We also verified 
that its energy is the same, within estimated numerical 
errors, for values of A lying within the chosen range. 
The binding energy of the first excited state of trimer 
is thus found to be 2.270 mK. We present this result, 
together with results of other sophisticated calculations 
reported so far in Table IV. 



TABLE IV: Results for energy of the first excited state of 
■'He-trimer (in units of mK) obtained by different methods. 



Reference 


Energy 


Ref [11] 


-2.282 


Ref [7] 


-2.280 


Ref [8] 


-2.277 


Present Method 


-2.270 



Here we remark that the isospectral formalism is an ef- 
ficient tool to calculate near-zero energy states in shal- 
low potential which supports at least one bound state. 
In the supersymmetric isospectral formalism, one can in 
principle calculate the energy value as well as the wave 
function. In our earlier calculation of the resonances of 
halo nuclei, we used the wave function of the isospec- 
tral potential. We calculated the probability of the sys- 
tem to be trapped in the well-barrier combination which 
facilitates the accurate determination of resonance en- 
ergy [13, HH . We have observed that the resonance en- 
ergy is independent of A. However in the present ap- 
plication to calculate the bound state wave function in 
the isospectral potentials and to calculate other physical 
obscrvables we should take a different approach. From 
Fig. 1, we see that decreasing A gradually keeps the long- 
range part of the isospectral potentials almost unchanged 
whereas the short-range part has a drastic change. Nat- 
urally the wave functions of the isospectral family are 
now A dependent and obviously the calculated physical 
obscrvables [e.g. average radius of the system] will be A 
dependent. In the SUSY isospectral formalism [l^, it 
is possible to calculate the wave function of the original 
potential a;o(r) from the wave functions of the isospectral 
potential a)o(^: The ground state Coif, A) and the first 
excited state Ci('', A) of the isospectral potential uJo{r,X) 
corresponding to the energies Eq and Ei respectively are 
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related to the ground state (o{r) and the first excited 
state Ci {r) of tlie original potential (r) by [l^ 



Co(r-,A) = 



Co(0 
/o + A 



and 



Ci(r,A) = (i?i-i?o)Ci+CoW(Co,Ci) 



(19) 



(20) 



where W(Co, Ci) = CoCi ~ CiCo is the Wronskian and Co('') 
and Ci(^) 8.re well-behaved functions. We can calculate 
the ground state wave function (o{r, A) of the isospectral 
potential ljq {r, A) from the original gound state wave 
function Co(^) (which is known as we calculated it to 
construct the family of isospectral potentials) using 
Eq. (19). Alternatively we can obtain (o{r,X) directly 
by solving for the isospectral potential Ldo{r,X). Also 
we can obtain (r, A) by solving for the isospectral 
potential (i>o(r. A) with energy Ei. Thus Eq. (20) is a 
first order differential equation in (r) and we can solve 
it subject to proper boundary conditions Ci(0) = a-iid 
C((0) 7^ to obtain Ci('')- With these wave functions we 
can calculate any physical properties of the system and 
this time they will be independent of the A. However 
usually in the SUSY isospectral formalism, people are 
interested only in the bound state energy and in this 
paper we are also interested only in the energy of the 
elusive first excited state. 

Next we discuss in detail about the choice of A 
parameter. As in our earlier works jsil - lssj . we observe 
that the energy is independent of the parameter A. But 
making A too small, may create large numerical error in 
the numerically calculated wave function Ci(^j A), as the 
well becomes extremely narrow and very deep. Then 
the wave function changes very rapidly over a very small 
interval. Clearly its derivatives arc inaccurate, which 
in turn affects the accuracy of the wave function at the 
next mesh point. The errors accumulate. As we solve 
the isospectral potential numerically, large numerical 
error will result for very small A. This may mask the 
overall accuracy. Thus the accuracy of the results are 
very crucially dependent on the choice of A. There is no 
prescribed rule to choose A. It is in general chosen by 
trial. One can choose finer mesh interval within the LSW 
for a small A. But the cumulative error at successive 
mesh point, increases as the total number of mesh points 
increases. Using the error estimate of the integration 



procedure (we use Runga-Kutta algorithm) we optimize 
the values of A and mesh interval to minimize the error. 
In our earlier calculations [sil [s^] for 2^ state of ^He, 
the shallowness of the effective potential is removed for 
A = 10, whereas A = 0.1 was required for 2^ state of 
the same system to get the expected behavior of the 
isospectral potential. Thus the optimum choice of A 
depends on the choice of the system and its state. 

Now in principle the ground state energy can also be 
recalculated using the isospectral potentials, in which 
the ground state of the original potential now lies in 
the very deep well of the isospectrals. But, when we 
solve the isospectral potential numerically, large error 
may creep due to the extreme narrowness of the well, 
as discussed above. Thus when we use the isospectral 
formalism for the real many-body problems, ground 
state energy is more accurate when determined by the 
original potential and the excited states and resonance 
states are more accurate when we solve the isospectral 
potentials. 



IV. CONCLUSIONS 

In conclusion, we remark that the isospectral potential 
with a judiciously chosen value of A can be very useful 
for an accurate calculation of near-zero energy states in 
a shallow potential and also the resonance states of halo 
and highly unstable systems. The present application to 
the very weakly bound and highly elusive first excited 
state in ^Hc trimer demonstrates the novelty and 
practical utility of this technique. 
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